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Abstract 

We develop quenched chiral perturbation theory for baryons using the graded- 
symmetry formalism of Bernard and Golterman and calculate non-analytic con- 
tributions to the baryon masses coming from quenched chiral loops. The usual 
term proportional to m^ 2 is substantially altered due to the cancellation of di- 
agrams with internal quark loops. In addition, the rf "hairpin" vertex leads to 
a new correction, proportional to m l J 2 . We compare our results to numerical 
lattice data and use them to estimate the size of the quenching error in the octet 
baryon masses. 
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1 Introduction and Formalism 



There have been several recent calculations of non-analytic loop corrections, or chiral logs, 
for the masses and amplitudes of quenched pseudo-Goldstone bosons |], • The aim has been 
to study the effect of the quenched approximation (QA) as the chiral limit is approached. 
In this note, we extend this work to baryons, using the Lagrangian approach of Bernard and 
Golterman |]. 

The key observation of Ref. |1[ is that quenched QCD possesses a larger chiral symmetry 
than the full theory. Following Morel, the fermionic determinant is removed from the func- 
tional integral by adding three pseudo-fermionic (i.e. commuting, spin-1/2) ghost quarks u, 
d and s, degenerate in mass with the physical quarks and coupled to the gauge fields in the 
same way. The resulting Lagrangian is symmetric, in the chiral limit, under the graded chiral 
symmetry SU(3|3)i x SU(3|3)#, and this symmetry must be used to construct a low-energy 
chiral Lagrangian. 

We introduce baryons using the "heavy-quark" formalism employed by Jenkins and 
Manohar in ordinary chiral perturbation theory (ChPT) ||. Our field B v satisfies all the 
kinematic properties of the baryon field B v of ref. || — i.e. B v (x) = exp(imi3pv^x^)B(x), 
B v = \(1 +y)B v , and S%B V = —\B V) where is the 4-velocity, and S£, the covariant "spin" 
operator, satisfies v -S v = 0. Henceforth we drop the subscript v. 

We find that, in the QA, the field B must be a three index tensor, Bijk, as opposed to 
the usual 3x3 matrix B. In an ordinary chiral theory, B \— > UBU\ where U is introduced 
by defining £(x) = exp(iir(x) / f) = S 1 / 2 (x), such that £ \— > L^W = U£RJ. The latter two 
relations are "quenched" by replacing 7r with the graded field $, defined by eq. (3) of ref. 
[0]; then L, R e SU(3|3)x, ! fl. To construct the baryons, we introduce the "Quark" vector 
Q = (u,d,s,u,d,s) which transforms under SU(3|3)y (i.e. L = R = U) as Qi \— > U^Qj. The 
baryon field is then defined to transform as a direct product of Quark vectors: 



K?7 



QTQfQT - Q?' a QTQk e abc (c^u, (i; 



where C = i^lo is the charge conjugation matrix. Generalizing to chiral transformations 
U(L, R, $), we obtain 

BJ jk - (-lf^^^^U^Uu.B}^ , (2) 

where we use a standard notation for the grading factor. (Indices corresponding to physical 
quarks (1-3) are read as 1, and those corresponding to ghost quarks (4-6) are read as 0.) 
The relevant part of the quenched chiral Lagrangian including baryons is C®> = Cbg + 
where 



P 



str(^E^E f ) + 2fistr(M + )\ + a d^ d^ - mfol (3) 
is the Bernard-Golterman Lagrangian, and 



C m = i(Bv-VB) + 2a(BS^BA fl ) + 2(3(BS^A^B) + 2-f(BS^B)str(A^) (4) 
+ a M (BBM + ) + p M (BM + B) + a(BB)sti(M + ). 
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Here, str is the supertrace, and the additional fields are defined as 

$ = ^str$ = j=( v '-f}'), 

where M = diag(m u , m d , m s , m u , m d , m s ). Our convention is such that /„- = 93MeV. In 
the brackets ( ) denote group-invariant contractions. (Anti-)symmetry under j <-> k in B^ 
leaves only two unique couplings to fields such as Aj , and these are denoted by the ordering. 

In the QA the rj' (and if) remain light, and the propagator is to be computed only from 
the S terms in £bg- Terms quadratic in the singlet $o are treated as additional (so-called 
hairpin) vertices; the insertion of the one proportional to 5 = mQ/(A8ir 2 f 2 ) leads to enhanced 
chiral logs. 

One subtlety in doing calculations is that there is a several-to-one mapping between the 
Bijk and the independent fields. In the physical (quark) sector Bijk\ q , where indices run from 
1 to 3, we find 

Bijk\q = -^{ £ ijk'B^ +Eikk>B J k ). (5) 

The Clebsch for a meson loop i on baryon line b can be written as 

C b,i = T,^a D %4- ( 6 ) 

The i\) h are normalization factors which can be obtained from @. For example, for the 
proton, 

V = (^112,^121,^211) = Vs(M,-2). (7) 

The D b ' 1 are numbers computed from the "flavor" part of the Feynman rules and can be 
represented as a sum of quark-flow diagrams, as illustrated below. 



2 Mass Renormalization 

The important loop corrections to the baryon masses mj are shown in Fig. [T|. The axial 
vertex (•) is a sum of three terms, "a," and "7," corresponding to the couplings in 
Cb<s>, and there is a similar sum for the symmetry-breaking mass vertex, denoted by the solid 
triangle. The baryon propagator can be written as 
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where the flavor diagrams denote delta functions in the matched indices and carry Clebsch 
and grading factors which we have suppressed. The pion propagator is diagonal in the QiQj 
basis of states. In the neutral sector, this is different from ordinary ChPT, where instead 
the propagator is diagonal in the 7r°, rj, rf basis. This difference is a manifestation of how the 
graded symmetry "quenches" the v( and forces the introduction of the hairpin vertex. 
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Figure 1: Baryon mass renormalization in quenched ChPT. • is the axial vertex 
solid triangle is the mass vertex («m, Pm, 0"), and x is the hairpin vertex (m^j 



(a, /3, 7), the 
= <548tt 2 / 2 ). 



Consider the term from Fig. [l|a proportional to (3 2 . It can be written as 

-i&m b = ip 2 Yl Cb,iX{Mi), (8) 

i=n,K,ss 

where X(M) = M 3 /(32irf 2 ) comes from the Feynman integral, and the Clebsch factor can 
be computed in terms of quark flow diagrams, as described above. For instance, if b denotes 
the proton, then the factors which enter into Eq. (§) are of the form 
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As expected, the graded symmetry produces a cancellation of all "quark loop" diagrams 
(the first term above is given by | Z)z={i, 1,0,0} ( — where the grading arises from the pion 
propagator), but there remain non-zero contributions from diagrams with only "valence 
quark" lines. Thus the mass renormalization proportional to ~ m^ 2 is substantially 
altered from its value in ordinary ChPT. The quenched nucleon, for example, will not be 
renormalized by kaon loops. 

Next, consider the hairpin insertion of Fig. |I|b. This Feynman integral contains an addi- 
tional pion propagator which reduces the order of the finite part by a power of M 2 . Thus 
a new non- analytic term (one proportional to M n ~ m l J 2 ) appears in the quenched expan- 
sion for baryon mass renormalization. The flavor factor in this case receives non-vanishing 
contributions only from the a and (3 vertices; terms containing a 7 vertex always have an 
internal quark loop. For instance, 



P 2 



a/3 



Pi 



We have calculated the renormalization for the octet baryons up to 0(M%) for m u = 
irid 7^ m s . We first present the result for degenerate quarks. To facilitate the comparison 
with ordinary ChPT, we insert © into and compute the correspondence between our 
parameters a, P, «m, etc. and the canonical ones defined in [||]. We find a = 2(^D + F), 



p = (~D+F), a M = 4(|6 D +6 F ), p M = 2(-|6 D +6 F ) and a = 2(b +b D - 
is unique to the quenched theory. Our result for the baryon mass is 



-bp). The 7- vertex 



m N 
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3F) 2 5il4 + 2(b D -3b F )M 2 + 
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\{D - 3F)(2D + 3 7 ) + \a Q {D - 3Ff 
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(9) 



3 



where Ml = 2fim q . There is a logarithm from Fig. lc which is absorbed into the tree- level 
result through pion mass renormalization: 



Additional contributions of order M% lnM^ are found to cancel. The cto term comes from 
the kinetic hairpin insertion, which is diagrammatically analogous to Fig. lb. 

3 Applications and Discussion 

We wish to use our results to estimate the errors due to quenching. Baryon masses have the 
advantage that the leading non-analytic corrections are larger than the analytic ones by a 
power of M n . (For meson masses, the enhancement is only logarithmic.) This absence of 
logarithms leads to the further advantage that no unknown scale A x is introduced. Thus, 
if we knew the quenched constants a, [3 etc., we could use the difference between quenched 
and full non-analytic terms as an estimate of the effect of quenching. Of course, in practice 
we do not know the constants; we must either assume that they are similar to their values 
in the full theory or find quantities in which they cancel. 

Such estimates will only be useful if the quark masses used in simulations are small enough 
for ChPT to be reliable. We first examine this question using the high statistics results for 
rriN from Ref. 0]. In Fig. 2 we plot tun vs. Ml at (3 = 5.93 (scaled using a -1 = 1.63) 
and make a "fit," with zero degrees of freedom, to the functional form obtained from Eq. 
(|Sp. We compare this to the "prediction" of our quenched calculation, using the full theory 
parameters from the fit of Ref. @], together with the assumptions 7 = «o = 0. We find, in 
units of GeV, 



We conclude from this comparison that ChPT is doing a reasonable job. Most significant 
is that the negative curvature predicted by ChPT, due to the Ml term, is seen with the 
correct magnitude. The "hook" at small m q is due to the M n term, which is absent in the 
full theory. The rough agreement on the sign and magnitude of this term is fortuitous, since 
both could be altered by adding in corrections. 

Proceeding next to quantitative estimates of quenching errors in baryon masses, we ad- 
dress the following question: if we could calculate dimensionless ratios (e.g. m N / f n ) in the 
quenched and full theories, using the physical quark masses]^ by how much would these ratios 
differ? The use of dimensionless ratios allows us to avoid worrying about overall differences 
in scale. First, consider the ratio m^j 'f n . The key observation here is that kaon loops do 
not contribute in the QA: 



Ml = 2^m q [l-8hi{M 2 JK)]- 



(10) 



ChPT: m N = 0.97 — 0.5— M w + 3.4 Ml -1.5M£, 
fit: m N = 0.96 - 1.0 M n + 3.6 Ml - 2.0M^. 




[ 1 + pion terms + . . . ] , 



3 We actually consider the slightly unphysical case of m. 
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Figure 2: The nucleoli mass from the lattice [6]. The curve is from a fit to the form = 

m + a(Ml) 1 ' 2 + b(M 2 ) + c(M 2 ) 3 / 2 . 



whereas they do in the full theory: 

i 

[1 + c 2 M 2 K - c 3 M| + c 4 M|ln(M| r /Aj) + pion terms + ...]. 

The coefficients q can be obtained from the literature. The chiral expansion in the two 
theories is clearly quite different. To make a quantitative estimate, we ignore the pion terms 
(for physical quark masses, the kaon terms are dominant), we assume that (jn/f)^ ¥ > ~ 
(m/ fY®' , and we evaluate the c, using the fit of [[|. We find that the ratio is 16% too large 
in the quenched approximation. If we set the scale using f n , this corresponds to a 150 MeV 
overestimate of m^. The error could be smaller if the difference in the chiral loops were 
compensated by a change in fn/ f, but this could only occur at a single value of M#. 

This estimate is about as good as one could hope to make, because only the lowest order 
coefficients m and / appear in the quenched result. We need not make any assumptions 
about the higher order coefficients etc. 

We can now do a similar analysis for all the members of the baryon octet. Rather than 
computing the ratios to f n , let us instead consider ratios within the octet itself, i2y = rrii/mj. 
These quantities have the advantages that (a) they do not involve the unknown scale A x , 
and (b) if we take their ratio in the full and quenched theories, then the tree-level mass terms 
cancel.^ As above, we consider physical quark masses and therefore drop terms from pion 
loops, which are numerically small. This leaves contributions from loops containing kaons 

4 They cancel assuming that we ultimately equate the quenched theory parameters with their full theory 
counterparts, as we have done above, and which we must do to obtain quantitative estimates. Henceforth 
this assumption will be implied, and we drop the (Q) and (F) superscripts. 




5 



and ss-mesons, the latter type being replaced by the rj in the full theory. In the QA, we find 
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where we have excluded both the O(ikf^) terms, which cancel in the final ratios, and the 
M\ terms proportional to ao, which vanish under our usual assumption that a§ ~ 0. In the 
full theory, the corresponding results are || || 
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The ratios = R^j/R^p determine the relative quenching error in the octet mass splittings. 
Given the assumptions stated above and taking the full-theory parameters from the first fit 
of Ref. [||, we obtain the following estimates: 

r^N = 1 + 0.19(5/0.2) + 0.13 = 1.31 [1.22] for S = 0.2 [0.1] , 
r BN = 1 - 0.46(5/0.2) + 0.43 = 0.97 [1.20] for 5 = 0.2 [0.1] , 
r AN = 1 - 0.39(5/0.2) + 0.26 = 0.87 [1.06] for S = 0.2 [0.1] . 

Clearly we need to know 5 to draw quantitative conclusions. The expectation using con- 
tinuum parameters is 5 = 0.2 |2[j, while the direct evidence suggests a smaller number 0. 
Nevertheless, the results indicate that the pattern of splittings within the octet might differ 
by 20-30% in the quenched approximation from that in the full theory. 

How do these estimates compare to numerical results? We can only make the compar- 
ison for m^v/ivr, since the most accurate baryon mass calculations have only been done for 
degenerate quarks. Ref. || || finds that this ratio is about 15% too large, in (overly) good 
agreement with our result. They find, however, that the p and decuplet baryon masses are 
also about 15% too large when compared to f n , and they suggest a picture in which the 
baryon masses are in the correct ratio, but that f n is too small. This hypothesis is in conflict 
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with our estimates for the remainder of the octet masses. It will be very interesting to see 
how simulations of quenched baryons with non-degenerate quarks turn out. 

There is one caveat in the above comparison: Ref. || obtains their value for using a 
linear extrapolation in m q , whereas our calculation includes terms proportional to m l J 2 and 
m?J 2 (see Fig. §). In other words the quenched computation of Ref. || is providing a quantity 

slightly different from rn§ calculated in quenched CHPT. It should be straightforward to 
account for this difference, and we will present results in a future write-up ||. 

Finally, we note that a shortcoming in the calculation discussed above is that it does not 
include decuplet baryons, neither as external states nor in loops. In the continuum, it has 
been found that the chiral expansion is, in general, better behaved if one includes the decuplet 
J], [3| . This is likely to be even more true when trying to describe lattice results, because the 
octet-decuplet splitting in the QA is considerably smaller. It is also interesting to work out 
the chiral expansions for the decuplet masses, since they are computed in present simulations 
and allow further tests of our methodology. We have nearly completed these calculations 
and will present results in ||. One result of note is that we find that the quenched A can 
decay to p + ir. 

Acknowledgements 

We would like to thank Claude Bernard, Maarten Golterman, Paul Mackenzie, Don Wein- 
garten, and Larry Yaffe for useful conversations. This work is supported in part by the DOE 
through grant DE-FG06-91ER40614, and by an Alfred P. Sloan Fellowship. 

References 

[1] Claude Bernard and Maarten Golterman, Phys. Rev. D 46, 853 (1992). 
[2] Stephen R. Sharpe, Phys. Rev. D 41, 3233 (1990). 

[3] Elizabeth Jenkins and Aneesh V. Manohar, Phys. Lett. B 255, 558 (1991). 
[4] V. Bernard, N. Kaiser and U. Meifiner, Z. Phys. C60, 111 (1993). 

[5] F. Butler, H. Chen, A. Vaccarino, J. Sexton and D. Weingarten, Phys. Rev. Lett. 70, 
2849 (1993). 

[6] Elizabeth Jenkins, Nucl. Phys. B 368, 190 (1992). 

[7] Y. Kuramashi, M. Fukugita, H. Mino, M. Okawa and A. Ukawa, these proceedings. 

[8] F. Butler, H. Chen, J. Sexton, A. Vaccarino, D. Weingarten, Meson Decay Constants 
from the Valence Approximation to Lattice QCD, IBM-HET-93-3, (Oct 1993), |hep- 
| lat/9310009| 

[9] J. Labrenz and S. Sharpe, in preparation. 



7 



